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Abstract 

The top quark mass was measured in the channels tt lepton-f-jets and tt dilepton 
(lepton=e, //) based on ATLAS data recorded in 2011. The data were taken at the LHC with a 
proton-proton centre-of-mass energy of = 7 TeV and correspond to an integrated lumin¬ 
osity of 4.6 fb”h The tt lepton-i-jets analysis uses a three-dimensional template technique 
which determines the top quark mass together with a global jet energy scale factor (JSF), 
and a relative Z^-to-light-jet energy scale factor (bJSF), where the terms Z>-jets and light-jets 
refer to jets originating from ^-quarks and u, d, c, 5-quarks or gluons, respectively. The ana¬ 
lysis of the tt dilepton channel exploits a one-dimensional template method using the nift, 
observable, defined as the average invariant mass of the two lepton-i-Z?-jet pairs in each event. 
The top quark mass is measured to be 172.33 + 0.75(stat -i- JSF -i- bJSF) + 1.02(syst) GeV, 
and 173.79 + 0.54(stat) + 1.30(syst) GeV in the tt lepton-i-jets and tt —> dilepton chan¬ 
nels, respectively. The combination of the two results yields mtop = 172.99 + 0.48(stat) + 
0.78(syst) GeV, with a total uncertainty of 0.91 GeV. 
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1. Introduction 


The mass of the top quark (mtop) is an important parameter of the Standard Model (SM) of particle 
physics. Precise measurements of mtop provide critical inputs to fits of global electroweak parameters [1- 
3] that help assess the internal consistency of the SM. In addition, the value of mtop affects the stability of 
the SM Higgs potential, which has cosmological implications [4-6]. 

Many measurements of mtop were performed by the CDF and DO collaborations based on Tevatron 
proton-antiproton collision data corresponding to integrated luminosities of up to 9.7 fb“^ A selection 
of these measurements was used in the recent Tevatron nitop combination resulting in mtop = 174.34 ± 
0.37 (stat) ± 0.52 (syst) GeV = 174.34 ± 0.64 GeV [7]. Since 2010, measurements of mtop from the LHC 
by the ATLAS and CMS collaborations have become available. They are based on proton-proton {pp) 
collisions at a centre-of-mass energy of = 7 TeV, recorded during 2010 and 2011 for integrated lumin¬ 
osities of up to 4.9 fb“^ [8-13]. The corresponding LHC combination, based on = 7 TeV data and in¬ 
cluding preliminary results, yields nitop = 173.29±0.23 (stat)±0.92(syst) GeV = 173.29±0.95 GeV [14]. 
Using the same LHC input measurements and a selection of the mtop results from the Tevatron experi¬ 
ments, the first Tevatron-tLHC mtop combination results in mtop = 173.34 + 0.27 (stat) ± 0.71 (syst) GeV, 
with a total uncertainty of 0.76 GeV [15]. Recently, improved individual measurements with a total un¬ 
certainty compatible with that achieved in the Tevatron-i-LHC nttop combination have become available; 
the most precise single measurement is obtained by the DO Collaboration using tt lepton-i-jets events 
and yields mtop = 174.98 ± 0.76 GeV [16]. 

This article presents a measurement of mtop using events with one or two isolated charged leptons (elec¬ 
trons or muons) in the final state (the tt —> lepton-i-jets and tt —> dilepton decay channels), in 4.6 fb“^ 
of pp collision data collected by the ATLAS detector at a centre-of-mass energy of = 7 TeV dur¬ 
ing 2011. It supersedes Ref. [8], where, using a two-dimensional fit to reconstructed observables in the 
tt lepton-i-jets channel, mtop was determined together with a global jet energy scale factor. The use 
of this scale factor allows the uncertainty on mtop stemming from imperfect knowledge of the jet energy 
scale (JES) to be considerably reduced, albeit at the cost of an additional statistical uncertainty compon¬ 
ent. The single largest systematic uncertainty on m^p in Ref. [8] was due to the relative (7-to-light-jet en¬ 
ergy scale (bJES) uncertainty, where the terms (?-jets and light-jets refer to jets originating from ( 7 -quarks 
and u, d, c, ^-quarks or gluons, respectively. To reduce this uncertainty in the present analysis, a three- 
dimensional template fit is used for the first time in the tt —> lepton-i-jets channel, again replacing the 
corresponding uncertainty by a statistical uncertainty and a reduced systematic uncertainty. This concept 
will be even more advantageous with increasing data luminosity. In addition, for the combination of the 
measurements of mtop in the two decay channels an in-depth investigation of the correlation of the two 
estimators for all components of the sources of systematic uncertainty is made. This leads to a much 
smaller total correlation of the two measurements than what is typically assigned, such that their combin¬ 
ation yields a very significant improvement in the total uncertainty on mtop. To retain this low correlation, 
the jet energy scale factors measured in the tt lepton-f-jets channel have not been propagated to the 
tt —7 dilepton channel. 

In the tt —7 lepton-i-jets channel, one W boson from the top or antitop quark decays directly or via an 
intermediate t decay into an electron or muon and at least one neutrino, while the other W boson decays 
into a quark-antiquark pair. The tt decay channels with electrons and muons are combined and referred 
to as the lepton-i-jets (or as a shorthand (?-i-jets) final state. The tt dilepton channel corresponds to 
the case where both W bosons from the top and antitop quarks decay leptonically, directly or via an 
intermediate t decay, into an electron or muon and at least one neutrino. The tt decay channels ee, ep,pp 
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are combined and referred to as the dilepton final state. For both the ^+jets and dilepton final stales, Ihe 
measuremenls are based on Ihe lemplale melhod [17]. In Ihis lechnique, Monle Carlo (MC) simulated 
disfribulions are conslrucfed for a chosen quanfify sensilive lo fhe physics paramefer under sludy, using 
a number of discrete values of lhaf paramefer. These femplafes are filled lo analylical functions lhal 
interpolate belween differenl inpul values of Ihe physics parameter, fixing all olher parameters of Ihe 
functions. In Ihe final step a likelihood fil lo Ihe observed dislribulion in data is used to obtain the value 
for the physics parameter that best describes the data. In this procedure the top quark mass determined 
from data corresponds to the mass definition used in the MC simulation. It is expected that the difference 
belween Ihis mass definition and Ihe pole mass is of order 1 GeV [18-21]. 

In Ihe T’+jels channel, evenls are reconslrucled using a kinematic fil lhal assumes a tt topology. A three- 
dimensional template method is used, where nitop is determined simultaneously with a light-jet energy 
scale factor (JSF), exploiting the information from the hadronic W decays, and a separate fj-to-light-Jet 
energy scale factor (bJSF). The JSF and bJSF account for residual differences of dala and simulation in Ihe 
lighl-jel and in Ihe relative fi-lo-lighl-jel energy scale, respectively, Ihereby mitigating Ihe corresponding 
systematic uncerlainlies on mtop. The analysis in Ihe dileplon channel is based on a one-dimensional 
template melhod, where Ihe templates are conslrucfed for Ihe mii, observable, defined as Ihe per-evenl 
average invarianl mass of Ihe Iwo leplon-fj-Jel systems from Ihe decay of Ihe lop quarks. Due lo Ihe 
underconslrained kinematics associated wilh Ihe dileplon final slate, no in silu conslrainl of Ihe jel energy 
scales is performed. 

This article is organised as follows: after a shorl description of Ihe ATLAS defector in Seel. 2, Ihe dala 
and MC simulation samples are discussed in Seel. 3. Delails of Ihe evenl selection and reconslruclion are 
given in Seel. 4. The template fils are explained in Seel. 5. The measuremenl of mtop in Ihe Iwo final states 
is given in Sect. 6 , and the evaluation of the associated systematic uncertainties are discussed in Sect. 7. 
The results of the combination of the mtop measurements from the individual analyses are reported in 
Sect. 8 . Finally, the summary and conclusions are given in Sect. 9. 


2. The ATLAS detector 


The ATLAS detector [22] covers nearly the entire solid angle around the collision point. * It consists of 
an inner tracking detector surrounded by a thin superconducting solenoid, electromagnetic and hadronic 
calorimeters, and a muon spectrometer incorporating three large superconducting toroid magnets. The 
inner-detector system (ID) is immersed in a 2 T axial magnetic field and provides charged-particle track¬ 
ing in the range \q\ < 2.5. The high-granularity silicon pixel detector covers the interaction region and 
typically provides three measurements per track, the first energy deposit being normally in the innermost 
layer. It is followed by the silicon microstrip tracker designed to provide four two-dimensional meas¬ 
urement points per track. These silicon detectors are complemented by the transition radiation tracker, 
which enables radially extended track reconstruction up to I 77 I = 2.0. The transition radiation tracker also 
provides electron identification information based on the fraction of energy deposits (typically 30 hits 
in total) above an energy threshold corresponding to transition radiation. The calorimeter system covers 


* ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the i/-axis points 
upwards. Cylindrical coordinates (r, 0) are used in the transverse plane, 0 being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 6 as 77 = - In tan(0/2). Angular distances are defined as A/? = 
V(A^)2 -H (A^)2. 
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the pseudorapidity range |? 7 l < 4.9. Within the region \r]\ < 3.2, electromagnetic calorimetry is provided 
by barrel and endcap high-granularity lead/liquid-argon (LAr) electromagnetic calorimeters, with an ad¬ 
ditional thin LAr presampler covering \q\ < 1.8, to correct for energy loss in material upstream of the 
calorimeters. Hadronic calorimetry is provided by the steel/scintillator-tile calorimeter, segmented into 
three barrel structures within \r]\ < 1.7, and two copper/LAr hadronic endcap calorimeters. The solid 
angle coverage is completed with forward copper/LAr and tungsten/LAr calorimeter modules optimised 
for electromagnetic and hadronic measurements respectively. The muon spectrometer (MS) comprises 
separate trigger and high-precision tracking chambers measuring the deflection of muons in the magnetic 
field generated by the toroids. The precision chamber system covers the region I 77 I < 2.7 with three layers 
of monitored drift tubes, complemented by cathode strip chambers in the forward region. The muon trig¬ 
ger system covers the range I 77 I < 2.4 with resistive plate chambers in the barrel, and thin gap chambers 
in the endcap regions. A three-level trigger system is used to select interesting events [23]. The Level-1 
trigger is implemented in hardware and uses a subset of detector information to reduce the event rate to at 
most 75 kHz. This is followed by two software-based trigger levels which together reduce the event rate 
to about 300 Hz. 


3. Data and Monte Carlo samples 

For the measurements described in this document, data from LHC pp collisions a.t ^[s - 1 TeV are 
used. They correspond to an integrated luminosity of 4.6 fb“' with an uncertainty of 1.8% [24], and were 
recorded during 2011 during stable beam conditions and with all relevant ATLAS sub-detector systems 
operational. 

MC simulations are used to model tt and single top quark processes as well as some of the background 
contributions. Top quark pair and single top quark production (in the s- and Wt-channels) are sim¬ 
ulated using the next-to-leading-order (NLO) MC program PowHEO-hvq (patch4) [25] with the NLO 
CTIO [26] parton distribution functions (PDFs). Parton showering, hadronisation and the underlying 
event are modelled using the Pythia (v6.425) [27] program with the Perugia 2011C (P2011C) MC para¬ 
meter set (tune) [28] and the corresponding CTEQ6L1 PDFs [29]. The AcerMC (v3.8) generator [30] 
interfaced with Pythia (v6.425) is used for the simulation of the single top quark t-channel process. The 
AcerMC and Pythia programs are used with the CTEQ6L1 PDEs and the corresponding P2011C tune. 

Eor the construction of signal templates, the tt and single top quark production samples are generated 
for different assumed values of m^op, namely 167.5,170,172.5,175,177.5 GeV. The tfMC samples are 
normalised to the predicted tt cross section for each mtop value. The tt cross section for pp collisions 
at = 7 TeV is cTtt = 177];}° pb for nitop = 172.5 GeV. It was calculated at next-to-next-to-leading- 
order (NNEO) in QCD including resummation of next-to-next-to-leading-logarithmic (NNEE) soft gluon 
terms with Top-i-i-2.0 [31-36]. The PDE-i-a^ uncertainties on the cross section were calculated using the 
PDE4EHC prescription [37] with the MSTW2008 68 % CE NNEO [38,39], CTIO NNEO [26,40] and 
NNPDE2.3 5f EEN [41] PDEs, and added in quadrature to the factorisation and renormalisation scale 
uncertainty. The NNEO-rNNEE value, as implemented in Hathor 1.5 [42], is about 3% larger than the 
plain NNEO prediction. The single top quark production cross sections are normalised to the approximate 
NNEO prediction values. Eor example, for nitop = 172.5 GeV, these are 64.6];2o pb [43], 4.6 ±0.2 pb [44] 
and 15.7 ± 1.1 pb [45] for the L, s- and Wt- production channels respectively. 

The production of VF or Z bosons in association with jets is simulated using the Alpgen (v2.13) gener¬ 
ator [46] interfaced to the Herwig (v6.520) [47,48] and Jimmy (v4.31) [49] packages. The CTEQ6E1 
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PDFs and the corresponding AUET2 tune [50] are used for the matrix element and parton shower set¬ 
tings. The VT-i-jets events containing heavy-flavour quarks {Wbb+]e.t?,, fPcc-i-jets, and VTc-i-jets) are gen¬ 
erated separately using leading-order matrix elements with massive b- and c-quarks. An overlap-removal 
procedure is used to avoid double counting of heavy-flavour quarks between the matrix element and 
the parton shower evolution. Diboson production processes {WW, VTZ and ZZ) are produced using the 
Herwig generator with the AUET2 tune. 

Multiple pp interactions generated with Pythia (v6.425) using the AMBT2B tune [51] are added to all 
MC samples. These simulated events are re-weighted such that the distribution of the number of interac¬ 
tions per bunch crossing (pile-up) in the simulated samples matches that in the data. The average number 
of interactions per bunch crossing for the data set considered is 8.7. The samples are processed through 
a simulation of the ATEAS detector [52] based on GEANT4 [53] and through the same reconstruction 
software as the data. 


4. Event selection and reconstruction 

4.1. Object selection 

In this analysis tt events with one or two isolated charged leptons in the final states are selected. The event 
selection for both final states is based on the following reconstructed objects in the detector: electron and 
muon candidates, jets and missing transverse momentum (E™®®). 

An electron candidate is defined as an energy deposit in the electromagnetic calorimeter with an asso¬ 
ciated well-reconstructed track [54]. Electron candidates are required to have transverse energy Ej > 
25 GeV and |?7ciusterl < 2.47, where Tyduster is the pseudorapidity of the electromagnetic cluster asso¬ 
ciated with the electron. Candidates in the transition region between the barrel and endcap calorimeter 
(1.37 < |?7ciusterl < 1-52) are excluded. Muon candidates are reconstructed from track segments in different 
layers of the MS [55]. These segments are combined starting from the outermost layer, with a procedure 
that takes effects of detector material into account, and matched with tracks found in the ID. The final 
candidates are refitted using the complete track information, and are required to satisfy pj > 20 GeV and 
I 77 I < 2.5. Isolation criteria, which restrict the amount of energy deposited near the lepton candidates, are 
applied to both the electrons and muons to reduce the backgrounds from heavy-flavour decays inside jets 
or photon conversions, and the background from hadrons mimicking lepton signatures, in the following 
referred to as non-prompt and fake-lepton background (NP/fake-lepton background). Eor electrons, the 
energy not associated with the electron cluster and contained in a cone of AR = 0.2 around the electron 
must not exceed an 77 -dependent threshold ranging from 1.25 to 3.7 GeV. Similarly, the total transverse 
momentum of the tracks contained in a cone of AR = 0.3 must not exceed a threshold ranging from 
1.00 to 1.35 GeV, depending on the electron candidate px and 77 . Eor muons, the sum of track transverse 
momenta in a cone of AR = 0.3 around the muon is required to be less than 2.5 GeV, and the total energy 
deposited in a cone of AR = 0.2 around the muon is required to be less than 4 GeV. The longitudinal 
impact parameter of each charged lepton along the beam axis is required to be within 2 mm of the recon¬ 
structed primary vertex, defined as fhe verfex wifh fhe highesf ZtrkPxtrk’ among all candidafes wifh af 
leasf five associated fracks wifh px.trk > 0.4 GeV. 

Jets are reconstructed with the anti-kt algorithm [56] using a radius parameter of = 0.4, starting from 
energy clusters of adjacent calorimeter cells called topological clusters [57]. These jets are calibrated first 
by correcting the jet energy using the scale established for electromagnetic objects (EM scale). They are 
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further corrected to the hadronic energy scale using calibration factors that depend on the jet energy and 77, 
obtained from simulation. Finally, a residual in situ calibration derived from both data and MC simulation 
is applied [58]. Jet quality criteria are applied to identify and reject jets reconstructed from energy deposits 
in the calorimeters originating from particles not emerging from the bunch crossing under study [59]. To 
suppress the contribution from low-px jets originating from pile-up interactions, tracks associated with 
the jet and emerging from the primary vertex are required to account for at least 75% of the scalar sum of 
the pt of all tracks associated with the jet. Jets with no associated tracks are also accepted. 

Muons reconstructed within a AR = 0.4 cone around a jet satisfying pj > 25 GeV are removed to 
reduce the contamination caused by muons from hadron decays within jets. Subsequently, jets within 
a AR = 0.2 cone around an electron candidate are removed to avoid double counting, which can occur 
because electron clusters are usually also reconstructed as jets. After this jet overlap removal, electrons 
are rejected if their distance to the closest jet is smaller than AR = 0.4. 

The reconstruction of is based on the vector sum of calorimeter energy deposits projected onto 
the transverse plane. The is reconstructed from topological clusters, calibrated at the EM scale 
and corrected according to the energy scale of the corresponding identified physics objects. Contribu¬ 
tions from muons are included by using their momentum as measured by the inner detector and muon 
spectrometer [60]. 

The reconstruction of top quark pair events is facilitated by the ability to tag jets originating from 
Z7-quarks. For this purpose the neural-network-based MVl algorithm is applied [61,62]. In the follow¬ 
ing, irrespective of their origin, jets tagged by this algorithm are called Z^-tagged jets, whereas those not 
tagged are called untagged jets. Similarly, whether they are tagged or not, jets originating from Z^-quarks 
and from u, d, c, ^-quarks or gluons are called Z^-jets and light-jets, respectively. The MV 1 algorithm re¬ 
lies on track impact parameters and the properties of reconstructed secondary vertices such as the decay 
length significance. The chosen working point corresponds to a f^-tagging efficiency of 75% for fj-jets in 
simulated ti events and a light-jet (c-quark jet) rejection factor of about 60 (4). To match the f7-tagging 
performance in the data, pt- and 77 -dependent scale factors are applied to MC jets depending on their 
original flavour. The scale facfors are obfained from dijef [62] and tt dilepfon evenfs. The 77-based 
calibration is obfained using fhe mefhodology described in Ref. [63], applied fo fhe 7 TeV dafa. The 
scale facfors are calculafed per jef and finally multiplied fo obfain an evenf weighf for any reconsfrucfed 
disfribufion. 


4.2. Event selection 

The tt lepfon-i-jefs signal is characferised by an isolafed charged lepfon wifh relafively high pj, E™®® 
arising from fhe neufrino from fhe lepfonic W boson decay, fwo f7-jefs and fwo lighf-jefs from fhe hadronic 
W boson decay. The main confribufions fo fhe background sfem from IT-i-jefs producfion and from fhe 
NP/fake-lepfon background. The normalisation of fhe JT-i-jefs background is esfimafed from dafa, based 
on fhe charge-asymmefry mefhod [64], and fhe shape is obfained from simulafion. For fhe NP/fake- 
lepfon background, bofh fhe shape of fhe disfribufions and fhe normalisafion are esfimafed from dafa by 
weighfing each selecfed evenf by fhe probabilify of confaining a NP/fake lepfon. This confribufion in bofh 
fhe election and fhe muon channel is esfimafed using a dafa-driven mafrix mefhod based on selecting 
fwo cafegories of evenfs, using loose and fighf lepfon selecfion requiremenfs [65]. The confribufions 
from single fop quark, Z-tjefs, and diboson producfion are faken from simulafion, normalised fo fhe besf 
available fheorefical cross secfions. 
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The tt —> dilepton events are characterised by the presence of two isolated and oppositely charged leptons 
with relatively high pT^ arising from the neutrinos from the leptonic W boson decays, and two 
Zj-jets. Background processes with two charged leptons from W- or Z decays in the final state, which 
are similar to the tt —> dilepton events, are dominated by single top quark production in the ITt-channel. 
Additional contributions come from Z+jets processes and diboson production with additional jets. In the 
analysis, these contributions are estimated directly from the MC simulation normalised to the relevant 
cross sections. Events may also be wrongly reconstructed as tt —> dilepton events due to the presence of 
NP/fake leptons together with fj-tagged jets and As for the tt lepton+jets channel, the NP/fake- 
lepton background is estimated using a data-driven matrix method [65]. 

The selection of tt event candidates consists of a series of requirements on the general event quality and the 
reconstructed objects designed to select events consistent with the above signal topologies. To suppress 
non-collision background, events are required to have at least one good primary vertex. It is required that 
the appropriate single-electron or single-muon trigger has fired; the trigger thresholds are 20 or 22 GeV 
(depending on the data-taking period) for the electrons and 18 GeV for muons. Candidate events in 
the ^-i-jets final stale are required to have exactly one reconstructed charged lepton with Et > 25 GeV 
for electrons, and px > 20 GeV for muons, matching the corresponding trigger object. Exactly two 
oppositely charged leptons, with at least one matching a trigger object, are required in the dilepton final 
state. In the p-i-jets channel, E™'*® > 20 GeV and E™*® -i- rri^ > 60 GeV are required.^ In the e-f-jets 
channel more stringent selections on and (E™®® > 30 GeV and > 30 GeV) are imposed 
due to the higher level of NP/fake-lepton background. Eor the ee and fin channels, in the dilepton final 
stale, E™** > 60 GeV is required. In addition, the invariant mass of the same-flavour charged-leplon 
pair, ma {tt = ee,iufi), is required to exceed 15 GeV, to reduce background from low-mass resonances 
decaying into charged lepton-antilepton pairs and Drell-Yan production. Similarly, to reduce the Z-rjets 
background, values of compatible with the Z boson mass are vetoed by requiring \ma - 91 GeV| > 
10 GeV. In the efi channel Hj > 130 GeV is required, where Hj is the scalar sum of the pj of the two 
selected charged leptons and the jets. Einally, the event is required to have at least four jets (or at least 
two jets for the tt dilepton channel) with pj > 25 GeV and \t]\ < 2.5. At least one of these jets must be 
fj-tagged for the tt lepton-f-jets analysis. In the dilepton final state, events are accepted if they contain 
exactly one or two fj-tagged jets. 

These requirements select 61786 and 6661 data events in the ti lepton-i-jets and tt dilepton channels, 
with expected background fractions of 22% and 2%, respectively. Due to their inherent mtop sensitivity, 
here and in the following, the single top quark processes are accounted for as signal in both analyses, and 
not included in the quoted background fractions. 


4.3. Event reconstruction 

After the event selection described in the previous section, the events are further reconstructed according 
to the decay topology of interest, and are subject to additional requirements. 


^ Here is the W boson transverse mass, defined as where the measured vector provides 

the neutrino (v) information. 




4.3.1. Kinematic reconstruction of the lepton+jets final state 


A kinematic likelihood fit [8, 66] is used to fully reconstruct the tt —> lepton+jets kinematics. The al¬ 
gorithm relates the measured kinematics of the reconstructed objects to the leading-order representation 
of the tt system decay. The event likelihood is constructed as the product of Breit-Wigner (BW) dis¬ 
tributions and transfer functions (TF). The W boson BW line-shape functions use the world combined 
values of the W boson mass and decay width from Ref. [3]. A common mass parameter, is used for 
the BW distributions describing the leptonically and hadronically decaying top quarks, and this is fitted 
event-by-event. The top quark width varies with and it is calculated according to the SM predic¬ 
tion [3]. The TF are derived from the Powheg +Pythia tt signal MC simulation sample at an input mass 
of mtop = 172.5 GeV. They represent the experimental resolutions in terms of the probability that the 
observed energy at reconstruction level is produced by a given parton-level object for the leading-order 
decay topology. 

The input objects to the likelihood are: the reconstructed charged lepton, the missing transverse mo¬ 
mentum and four jets. For the sample with one b-tagged jet these are the f^-tagged jet and the three un¬ 
tagged jets with the highest pj. For the sample with at least two fj-tagged jets these are the two highest-px 
fj-tagged jets, and the two highest-px remaining jets. The x- and ^-components of the missing transverse 
momentum are used as starting values for the neutrino transverse momentum components, with its longit¬ 
udinal component {pv,z) as a free parameter in the kinematic likelihood fit. Its starting value is computed 
from the VT —> fv mass constraint. If there are no real solutions for pv,^ a starting value of zero is used. If 
there are two real solutions, the one giving the largest likelihood value is taken. 

Maximising the event-by-event likelihood as a function of establishes the best assignment of recon¬ 
structed jets to partons from the tt —> lepton+jets decay. The maximisation is performed by testing all 
possible permutations, assigning jets to partons. The likelihood is extended by including the probability 
for a jet to be f^-tagged, given the parton from the top quark decay it is associated with, to construct an 
event probability. The f?-tagging efficiencies and rejection factors are used to favour permutations for 
which a fj-tagged jet is assigned to a fj-quark and penalise those where a f^-tagged jet is assigned to a light 
quark. The permutation of jets with the highest likelihood value is retained. 

The value of obtained from the kinematic likelihood fit is used as the observable primarily sensitive 
to the underlying mtop. The invariant mass of the hadronically decaying W boson is calculated 

from the assigned jets of the chosen permutation. Finally, an observable called designed to be 
sensitive to the relative fj-to-light-jet energy scale, is computed in the following way. For events with 
only one b-tagg&d jet, is defined as the ratio of the transverse momentum of the f?-tagged jet to the 
average transverse momentum of the two jets of the hadronic W boson decay. For events with two or more 
b-tagg&d jets, is defined as the scalar sum of the transverse momenta of the f^-tagged jets assigned 
to the leptonically and hadronically decaying top quarks divided by the scalar sum of the transverse 
momenta of the two jets associated with the hadronic W boson decay. The values of and are 
computed from the jet four-vectors as given by the jet reconstruction to keep the maximum sensitivity to 
changes of the jet energy scale for light-jets and f7-jets. 

In view of the template parameterisation described in Sect. 5 additional selection criteria are applied. 
Events in which a fj-tagged jet is assigned to the W decay by the likelihood fit are discarded. This is 
needed to prevent mixing effects between the information provided by the and distributions. 
The measured is required to be in the range 125 GeV to 225 GeV for events with one fj-tagged jet, 
and in the range 130 GeV to 220 GeV for events with at least two fj-tagged jets. In addition, is 
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required to be in the range 55 GeV to 110 GeV and finally, is required to be in the range 0.3 to 3.0. 
The fraction of data events which pass these requirements is 35%. Although removing a large fraction of 
data, these requirements remove events in the tails of the three distributions, which are typically poorly 
reconstructed with small likelihood values and do not contain significant information on mtop. In addition, 
the templates then have simpler shapes which are easier to model analytically with fewer parameters. 


4.3.2. Reconstruction of the dilepton final state 

In the tt dilepton channel the kinematics are under-constrained due to the presence of at least two 
undetected neutrinos. Consequently, instead of attempting a full reconstruction, the nitop-sensitive ob¬ 
servable mgh is defined based on the invariant mass of the two charged-lepton-i-fj-jet pairs. 

The preselected events contain two charged leptons, at least two jets, of which either exactly one or 
exactly two are f^-tagged. For events with exactly two f7-tagged jets the charged-lepton-i-f7-tagged jet pairs 
can be built directly. In the case of events with only one fj-tagged jet the missing second fj-jet is identified 
with the untagged jet carrying the highest MVl weight. For both classes of events, when using the two 
selected jets and the two charged leptons, there are two possible assignments for the jet-lepton pairs, each 
leading to two values for the corresponding pair invariant masses. The assignment resulting in the lowest 
average mass is retained, and this mass is taken as the estimator of the event. The measured 
is required to be in the range 30 GeV to 170 GeV. This extra selection retains 97% of the data candidate 
events. 


4.3.3. Event yields 

The numbers of events observed and expected after the above selections are reported in Table 1 for the 
^-i-jets and dilepton final states. The observed numbers of events are well described by the sum of the 
signal and background estimates within uncertainties. The latter are estimated as the sum in quadrature 
of the statistical uncertainty, the uncertainty on the f?-tagging efficiencies, a 1.8% uncertainty on the 
integrated luminosity [24], the uncertainties on the tt and single top quark theoretical cross sections, a 30% 
uncertainty on the VT-i-jets and Z-rjets normalisation, and finally a 50% uncertainty on the NP/fake-lepton 
background normalisation. The distribution of several kinematic variables in the data were inspected and 
found to be well described by the signal-plus-background prediction, within uncertainties. As examples. 
Fig. l(left) shows the distribution of the untagged and f?-tagged jets pj observed in the ^-i-jets final state. 
Similarly, the pj distributions for the charged leptons and f7-tagged jets in the dilepton final state are 
shown on the right of Fig. 1 . In all cases the data are compared with the MC predictions, assuming an 
input top quark mass of 172.5 GeV. 


5. Analysis method 

The observables exploited in the mtop analyses are: in the tt —> lepton-i-jets channel and 

in the tt dilepton channel. 

In the tt lepton-rjets channel, templates of are constructed as a function of the top quark mass 
used in the MC generation in the range 167.5-177.5 GeV, in steps of 2.5 GeV. In addition, for the central 
mass point, templates of mJJop” are constructed for an input value of the light-jet energy scale factor (JSF) 
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(a) tt —» lepton+jets: untagged jets from W boson decays 


(b) tt —> dilepton: charged leptons 
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(c) tt —> lepton+jets: ^-tagged jets 


(d) tt —» dilepton: Zt-tagged jets 


Figure 1: Distributions of the transverse momentum of the untagged and Z>-tagged jets in the tt —> lepton+jets 
analysis (a, c) and of the charged lepton and Z>-tagged jets pj in the tt —> dilepton analysis (b, d). The data are 
shown by the points and the signal-plus-background prediction by the solid histogram. The hatched area is the 
combined uncertainty on the prediction described in Sect. 4.3.3, and the rightmost bin contains the overflow if 
present. For each figure, the ratio of the data to the MC prediction is also presented. 
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^H-jets final state 

Process 

One Z7-tagged jet 

At least two Z7-tagged jets 

Sum 

tt signal 

9890 + 

630 

8210 + 


560 

18100 + 

1100 

Single top quark (signal) 

756 + 

41 

296 + 


19 

1052 + 

57 

IkH-jets (data) 

2250 + 

680 

153 + 


49 

2400 + 

730 

Z-Hjets 

284 + 

87 

18.5 + 


6.1 

303 + 

93 

WW/WZ/ZZ 

43.5 + 

2.3 

4.65 + 


0.48 

48.2 + 

2.6 

NP/fake leptons (data) 

700 + 

350 

80 + 


41 

780 + 

390 

Signal-Hbackground 

13920 + 

1000 

8760 + 


560 

22700 + 

1400 

Data 

12979 

8784 

21763 

Exp. Bkg. frac. 

0.25 + 

0.02 

0.03 + 


0.00 

0.16 + 

0.01 

Data/MC 

0.93 + 

0.07 

1.00 + 


0.07 

0.96 + 

0.06 

Dilepton final state 

Process 

One Z7-tagged jet 

Two Z7-tagged jets 

Sum 

tt signal 

2840 + 

180 

2950 + 


210 

5790 + 

360 

Single top quark (signal) 

181 + 

10 

82.5 + 


5.7 

264 + 

15 

Z-Hjets 

34 + 

11 

4.1 + 


1.5 

38 + 

12 

WW/WZ/ZZ 

7.01 + 

0.63 

0.61 + 


0.15 

7.62 + 

0.67 

NP/fake leptons (data) 

52 + 

28 

2.6 + 


8.4 

55 + 

30 

Signal-Hbackground 

3110 + 

180 

3040 + 


210 

6150 + 

360 

Data 

3227 

3249 

6476 

Exp. Bkg. frac. 

0.03 + 

0.00 

0.00 + 


0.00 

0.02 + 

0.00 

Data/MC 

1.04 + 

0.06 

1.07 + 


0.07 

1.05 + 

0.06 


Table 1; The observed numbers of events, according to the h-tagged jet multiplicity, in the ^+jets and dilepton 
final states in 4.6 fb ' of ^fs = 7 TeV data. In addition, the expected numbers of signal and background events 
corresponding to the integrated luminosity of the data are given. The predictions are quoted using two significant 
digits for their uncertainty. The MC estimates assume SM cross sections. The VT+jets and NP/fake-lepton back¬ 
ground contributions are estimated from data. The uncertainties for the estimates include the components detailed 
in Sect. 4.3.3. Values smaller than 0.005 are listed as 0.00. 


in the range 0.95-1.05 in steps of 2.5% and for an input value for the relative Z>-to-light-jet energy scale 
factor (bJSF) in the same range. Independent MC samples are used for the different nttop mass points, 
and from those samples templates with different values of JSF and bJSF are extracted by appropriately 
scaling the four-momentum of the jets in each sample. The input value for the JSF is applied to all jets, 
whilst the input value for the bJSF is applied to all f?-jets according to the information about the generated 
quark flavour. This scaling is performed afler fhe various correction sfeps of fhe jef calibration and before 
any even! selecfion. This resulfs in differenl evenfs enfering fhe final selecfion from one energy scale 
variation fo anofher. Similarly, femplafes of rrf^° are consfrucfed as a funcfion of an inpuf JSF combining 
fhe samples from all mtop mass poinfs. Finally, femplafes of are consfrucfed as a funcfion of rritop, 
and as a funcfion of an inpuf bJSF af fhe cenfral mass poinf. 

In fhe tt —> dilepfon channel, signal femplafes for are consfrucfed as a funcfion of fhe fop quark 
mass used in fhe MC generation in fhe range 167.5-177.5 GeV, using separate samples for each of fhe 
five mass poinfs. 

The dependencies of fhe and disfribufions on fhe underlying mtop used in fhe MC simulation 
are shown Fig. 2(a, b), for evenfs with at least (exactly) two Z^-tagged jets, for the tt lepton-i-jets 
(tt dilepton) channel. The and distributions shown in Fig. 2(c-f), exhibit sizeable sensitivity 
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to global shifts of the JSF and the bJSF. These effects introduce large systematic uncertainties on mtop 
originating from the uncertainties on the JES and bJES, unless additional information is exploited. As 
shown for the tt —> lepton+jets channel in Eig. 3(a, c, e), the distribution is sensitive to changes of 
the JSE, while preserving its shape under variations of the input mtop and bJSE. As originally proposed 
in Ref. [17], a simultaneous fit to and is used to mitigate the JES uncertainty. The 
distributions show substantial sensitivity to the bJSE, and some dependence on the assumed mtop in the 
simulation, Eig. 3(b, d, f). Complementing the information carried by the mf®™ and observables, 
y^reco jjj unbinned likelihood fit to the data to simultaneously determine mtop, JSE, and bJSE. 

The per-event correlations of any pair of observables (mj:®™, and are found to be smaller than 
0.15 and are neglected in this procedure. 

5.1. Templates and fits in the tt -* lepton+jets channel 

Signal templates are derived for the three observables for all mtop-dependent samples, consisting of the tt 

signal events, together with single top quark production events. The signal templates for the 

and distributions are fitted to the sum of a Gaussian function and a Eandau function for ml®™ and 

bq top 

7?^“°, and to a sum of two Gaussian functions for (Eigs. 2 and 3). Eor the background, the mf®™ 
distribution is fitted to a Eandau function, while both the and the distributions are fitted to the 
sum of two Gaussian functions. To exploit the different sensitivities to the underlying mtop, JSE and bJSE, 
all template fits are performed separately for events with one T^-tagged jet, and for events with at least two 
Tj-tagged jets. 

Erom individual fits to all signal templates listed above, it was verified that the parameters of the fitting 
functions depend linearly on the respective parameter mtop, JSE or bJSE. Consequently, this linearity 
is imposed when parametrising the fitting functions for the combined fit to all signal templates for the 
three observables. Eor the signal, the parameters of the fitting functions for mj®™ depend linearly on mtop, 
JSE and bJSE. The parameters of the fitting functions of m'^^° depend linearly on the JSE. Einally, the 
parameters of the fitting functions of depend linearly on the bJSE and on m^p. As shown in Eig. 3, 
the dependencies of nf^° on m^p and bJSE, and of on JSE are negligible. Eor the background, the 
parameter dependencies of the fitting functions are the same except that, by construction, they do not 
depend on mtop. 

Signal and background probability density functions P®*® and for the mj®p°, and distribu¬ 
tions are used in an unbinned likelihood fit to the data for all events, i = 1,... A. The likelihood function 
maximised is: 


N 



(mtop, JSE, bJSE, /bkg) = Y] I "^top, JSE, bJSE, Mg) 


i=\ 


xPw(m;^“’'|JSE,/bkg) 

xPfo,(P;f“’Vtop,bJSE,/bkg), 


^reco,i 


( 1 ) 


13 








Figure 2; Distributions of in the tt lepton+jets channel (left) and in the tt —> dilepton channel (right) 
and their template parameterisations for the signal, composed of simulated ff and single top quark production events. 
The expected sensitivities of and are shown for events with at least two (or exactly two) h-tagged jets. 
Figures (a, b) report the distributions for different values of the input Wtop (167.5, 172.5 and 177.5 GeV). Figures (c, 
d) and (e, f) show the and distribution for mtop =172.5 GeV, obtained with JSF or bJSF of 0.95, 1.00 and 
1.05, respectively. Each distribution is overlaid with the corresponding probability density function that is obtained 
from the combined fit to all signal templates for all abservables. 
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Figure 3: Distributions of (left) and (right) in the tt lepton+jets channel and their template paramet- 

erisations for the signal, composed of simulated tt and single top quark production events. The expected sensitivity 
of and are shown for events with at least two h-tagged jets. Figures (a, b) report the distributions for 
different values'^of the input mtop (167.5, 172.5 and 177.5 GeV). Figures (c, d) and (e, f) show the and 
distribution for mtop =172.5 GeV, obtained with JSF or bJSF of 0.95, 1.00 and 1.05, respectively. Each distribution 
is overlaid with the corresponding probability density function that is obtained from the combined fit to all signal 
templates for all abservables. 
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with: 


I /«top, JSF, bJSF, /bkg) - (1 - /bkg) • F,^^p«p°’'' I "^top, JSF, bJSF) + 

/bkg-F;5K:riJSF,bjsF), 

I JSF, /bkg) - (1 - /bkg) • | JSF) + 

/bkg-PS^Vr’'’|JSF), 

I /^top, bJSF, /bkg) - (1 - /bkg) • I + 

/bkg I bJSF) 

where the fraction of background events is denoted by /bkg- The parameters to be determined by the fit 
are nitop, JSF, bJSF and /bkg, where /bkg is determined separately for the tt —> lepton+jets data sets with 
exactly one or at least two Z^-tagged jets. 

Pseudo-experiments are used to verify the internal consistency of the fitting procedure and to obtain the 
expected statistical uncertainty corresponding to a data sample of 4.6 fb“^ For each choice of the input 
parameters, 500 pseudo-experiments are generated. To retain the correlation of the analysis observables, 
individual MC events drawn from the full simulated event samples are used, rather than sampling from 
the separate mfop”, and distributions. For all five parameters, good linearity is found between 
the input parameters used to perform the pseudo-experiments, and the results of the fits. Within their 
statistical uncertainties, the mean values and widths of the pull distributions are consistent with the ex¬ 
pectations of zero and one, respectively. This means the method is unbiased with appropriate statistical 
uncertainties. The expected statistical uncertainties on mtop including the statistical contributions from 
the simultaneous fit of the JSF and bJSF obtained from pseudo-experiments at an input top quark mass 
of mtop = 172.5 GeV, and for a luminosity of 4.6 fb“', are 1.50 ± 0.06 GeV and 0.89 + 0.01 GeV for 
the case of one ^-tagged jet and for the case of at least two Z^-tagged jets, respectively. The results cor¬ 
respond to the mean value and the standard deviation of the distribution of the statistical uncertainties of 
the fitted masses from the pseudo-experiments. The different expected statistical uncertainties on mtop for 
the samples with different numbers of fi-tagged jets, which are obtained from samples containing similar 
numbers of events (see Table 1), are mainly a consequence of the different resolution on mtop. 

5.2. Templates and fits in the tt dilepton channel 

The signal templates comprise both the tt and the single top quark production processes, and are 
fitted to the sum of a Gaussian function and a Landau function, while the background distribution is fitted 
to a Landau function. Similarly to the tt lepton-rjets channel, all template fits are performed separately 
for events with one fi-tagged jet, and for events with exactly two fi-tagged jets. In Fig. 2(b) the sensitivity 
of the observable to the input value of the top quark mass is shown for the events with exactly two 
Zj-tagged jets, by the superposition of the signal templates and their fits for three input mtop values. For 
the signal templates, the parameters of the fitting functions of depend linearly on mtop. 

Signal and background probability density functions for the estimator are built, and used in an 
unbinned likelihood fit to the data for all events, i = I,.. .N. The likelihood function maximised is: 

^shaTe™<^'”t°P’ /bkg) = n I "hop) + /bkg ' ’ (2) 

i=l 
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where, as for the tt lepton+jets case, Pj'® and are the signal and background probability density 
functions and f\±g is the fraction of background events in the selected data set. 

Using pseudo-experiments, also for this decay channel good linearity is found between the input top quark 
mass used to perform the pseudo-experiments, and the results of the fits. Within their statistical uncertain¬ 
ties, the mean values and widths of the pull distributions are consistent with the expectations of zero and 
one, respectively. The expected statistical uncertainties on nitop obtained from pseudo-experiments for an 
input top quark mass of mtop = 172.5 GeV, and for a luminosity of 4.6 fb“^, are 0.95 ± 0.04 GeV and 
0.65 + 0.02 GeV for events with exactly one or two Zj-tagged jets, respectively. As for the ^-rjets channel, 
the different expected statistical uncertainties on mtop for the samples with different numbers of f?-tagged 
jets, which are obtained from samples containing similar numbers of events (see Table 1), are mainly a 
consequence of the different resolution on mtop. 


5.3. Combined likelihood fit to the event samples 

The final results for both the f’-i-jets and dilepton final states are obtained combining at the likelihood level 
the events with one or more f7-tagged jets. The measured mtop is assumed to be the same in these two 
sub-samples per decay channel. Similarly, the JSF and the bJSF are taken to be the same for the samples 
of the tt —> lepton-i-jets analysis with different 6-tagged jet multiplicities. On the contrary, the background 
fractions for the two decay channels, and for the samples with different numbers of fj-tagged jets, are kept 
independent, corresponding to four individual parameters ■: fttg^ ’ 

The combined likelihood fit allows the statistical uncertainties on the fitted parameters to be reduced, 
while mitigating some systematic effects. The expected statistical precision on mtop, for an input top 
quark mass of mtop = 172.5 GeV, a luminosity of 4.6 fb“^, and in the combined one or more f^-tagged 
jets event sample, is 0.76 + 0.01 GeV and 0.54 ± 0.01 GeV for the tt lepton-rjets and tt dilepton 
analyses, respectively. 


6. Top quark mass measurements 


The results of the fits for the tt lepton-f-jets and tt dilepton analyses are: 


^+jets 


- 172.33 + 0.75 (stat -t JSF -t bJSF) GeV, 
JSF - 1.019 ±0.003 (stat), 
bJSF - 1.003 ±0.008 (stat), 

173.79 ± 0.54 (stat) GeV. 


m, 


dil 

top 


For the tt lepton±jets channel, the fitted background fractions amount to 18.4%±2.2% and2.4%±1.5% 
for one f?-tagged jet and the at least two f7-tagged jets samples respectively. The corresponding values for 
the tt dilepton analysis are 3.5% ±3.7% and 1.4% ±2.2% for one f7-tagged jet and the two f7-tagged jets 
samples respectively. All quoted uncertainties are statistical only. These fractions are consistent with the 
expectations given in Table 1. The correlation matrices for the fitted parameters in the tt —> lepton±jets 
and tt dilepton analyses are reported in Table 2. 

Figure 4 shows the and distributions in the data together with the corresponding 

fitted probability density functions for the background alone and for the sum of signal and background. 
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JSF bJSF XkT‘'’“ 

?+]efs 

"^top 

JSF 

bJSF 

rf’+jets,lfe 

7bkg 

J+jcts.lb 

7bks 

1.00 

-0.36 1.00 

-0.89 0.03 1.00 

-0.03 -0.01 0.06 1.00 

-0.06 -0.09 0.09 0.01 1.00 



dil ^clihlfj 

"^top Jhkg Jbkg 

"^op 

7bkg 

^il,2Z7 

Jhkg 

1.00 

0.07 1.00 

-0.14 -0.01 1.00 


Table 2: The correlations of the fitted parameters used in the likelihood maximisation of the tt lepton+jets 
analysis (left) and the tt —» dilepton analysis (right). 


The uncertainty bands are obtained by varying the three fitted parameters nitop, JSF, and bJSF within 
±lcr of their full uncertainties taking into account their correlation, while keeping the background frac¬ 
tions fixed. The individual systematic uncertainties and the correlations are discussed in Sect. 7 and 
Sect. 8, respectively. The band shown is the envelope of all probability density functions obtained from 
500 pseudo-experiments varying the parameters. Within this band, the data are well described by the 
fitted probability density function. 

For the tt —> lepton-i-jets analysis, the measured values of the three observables JSF, and bJSF), 

together with two-dimensional statistical uncertainty contours (± Icr), including the statistical components 
from the JSF and bJSF determination, are shown in Fig. 5(a-c). Correspondingly, the likelihood profile 
as a function of is reported in Fig. 5(d), for the sample with one Z^-tagged jet, the sample with two 
b-tagged jets and the combined tt dilepton result. These results demonstrate the good agreement 
between the parameter values measured in the samples with different fi-tagged jet multiplicities. 


7. Uncertainties affecting the mtop determination 

Several sources of systematic uncertainty are considered. Their effects on the ^-tjets and dilepton meas¬ 
urements are listed in Table 3, together with the result of the combination of the two channels discussed 
in Sect. 8. Each source of uncertainty considered is investigated, when possible, by varying the relev¬ 
ant quantities by ±Icr with respect to their default values. Using the changed parameters, 500 pseudo¬ 
experiments are performed using events drawn from the full simulated samples. The difference of the av¬ 
erage mtop computed from pseudo-experiments based on the standard MC sample, and the varied sample 
under consideration, both evaluated with the original template parameterisations, is used to determine the 
corresponding uncertainty. Unless stated otherwise, the systematic uncertainties arising from the different 
modelling sources are calculated as half of the difference of fhe resulfs of fhe upward and downward vari- 
afions. The sysfemafic uncertainties for fhe measured JSF and bJSF in fhe ti lepfon-i-jefs final sfafe are 
also esfimafed. Following Ref. [67], fhe acfual observed difference is quofed as fhe sysfemafic uncertainfy 
on fhe corresponding source, even if if is smaller fhan ifs associafed sfafistical precision. The latter is es- 
timafed faking info account the statistical correlation of the MC samples used in the comparison. The total 
uncertainty is calculated as the sum in quadrature of all individual contributions, i.e. neglecting possible 
correlations (small by construction). The estimation of the uncertainties for the individual contributions 
is described in the following. 
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Figure 4; The fitted distributions in the data, showing (a) m™, (b) (c) and (d) The fitted prob¬ 

ability density functions for the background alone and for signal-plus-background are also shown. The uncertainty 
bands indicate the total uncertainty on the signal-plus-background fit obtained from pseudo-experiments as ex¬ 
plained in the text. Figures (a-c) refer to the ff —> lepton-i-jets analysis, figure (d) to the ff —> dilepton analysis. 
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Figure 5; Likelihood contours showing the correlation determined in data of the measured to (a) the JSF 

and (b) the bJSF, and (c) the correlation of the two scales JSF and bJSF, within the tt —> lepton+jets analysis. 
Figures (a-c) show the results using the events with one ^-tagged jet only (grey ellipses), with at least two Z7-tagged 
jets (red ellipses) and hnally with all selected events, i.e. the ones with at least one /j-tagged jet (blue ellipses). 
The ellipses correspond to the +lcr (statistical) uncertainties, including the statistical components from the JSF and 
bJSF determination. While tracing the contours the additional parameters of the likelihood are hxed to their best 
ht values. Figure (d) reports the likelihood prohle as a function of for the sample with one /j-tagged jet, the 
sample with two Z7-tagged jets and the combined result. The colour coding is analogous to hgures (a-c). 
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7.1. Statistics and method calibration 


Statistical components due to the jet energy scale factors 

The statistical uncertainty quoted for the tt lepton+jets analysis is made up of three parts: a purely 
statistical component on mtop and the contributions stemming from the simultaneous determination of 
the JSF and bJSF. The former is obtained from a one-dimensional template method exploiting only the 
observable (fixing the values of the JSF and bJSF to the results of the three-dimensional analysis). 
The contribution to the statistical uncertainty on the fitted parameters due to the simultaneous fit of mtop 
and JSF, is estimated as the difference in quadrature of the statistical uncertainty of a two-dimensional 
and m^^°, fixing fhe value of bJSF) fif and fhe one-dimensional fif fo fhe dafa described above. 
Analogously, fhe confribufion of fhe sfafisfical uncerfainfy due fo fhe simulfaneous fif of bJSF fogefher 
wifh nitop and JSF, is defined as fhe difference in quadrafure of fhe sfafisfical uncerfainfies obfained in fhe 
fhree-dimensional and fhe fwo-dimensional (fixing bJSF) fits to the data. This separation allows a direct 
comparison of the sensitivity of the rriiop estimator for any analysis, irrespective of the number of observ¬ 
ables exploited by the fit. In addition, the sensitivity of the estimators for the global jet energy scales can 
be directly compared. These uncertainties can be treated as uncorrelated uncertainties in mtop combina¬ 
tions. Together with the systematic components of the residual jet energy scale uncertainty discussed in 
Sect. 7.4 below, they directly replace the uncertainty on mtop from the jet energy scale variations present 
without the in situ determination. 


Method calibration 

This uncertainty takes into account the effect of any bias introduced in the fit by the presence of correla¬ 
tions among the observables (neglected in the fit for the tt —> lepton-i-jets analysis) as well as the impact 
of the limited size of the MC samples (for both analyses). This leads to a systematic uncertainty in the 
template fit, which is reflected in the residual mass differences of fhe fiffed mass and fhe inpuf mass for a 
given MC sample. The largesf average difference observed in fhe pseudo-experimenfs carried ouf varying 
fhe underlying fop quark mass, fhe JSF and fhe bJSF wifh respecf fo fhe respecfive inpuf paramefer, is 
faken as fhe uncerfainfy from fhis source. 


7.2. tt modelling 

Signal Monte Carlo generator 

The sysfemafic uncertainty related to the choice of ti signal generator program is determined by com¬ 
paring the results of pseudo-experiments performed with either the MC@NLO [68,69] samples or the 
PowHEG samples, both generated with nitop = 172.5 GeV and using the Herwig program to perform the 
hadronisation. This choice is supported by the observation that these MC@NLO and Powheg samples 
exhibit very different jet multiplicities for the tt —> lepton-i-jets channel which bracket those observed 
in data [70]. The full difference of the results averaged over all pseudo experiments is quoted as the 
systematic uncertainty. 

The impact of changing the factorisation and renormalisation scales (jUf/r) in Powheg was also checked. 
The resulting mtop systematic uncertainties amount to 0.15 ± 0.07 GeV and 0.14 ± 0.05 GeV for the tt 
lepton-i-jets channel, and tt dilepton analysis respectively. Within the quoted statistical uncertainties, 
the jUp/R systematic uncertainties are consistent with those originating from the comparison of MC@NLO 
and Powheg, which are used here. 
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Hadronisation 

Signal samples for mtop = 172.5 GeV from the Powheg event generator are produced performing the 
parton showering and the hadronisation with either Pythia with the P2011C tune or Herwig and Jimmy 
with the ATLAS AUET2 tune [50]. The full difference of the results averaged over all pseudo experiments 
is quoted as the systematic uncertainty. 


Initial- and final-state QCD radiation 

Different amounts of initial- and final-state QCD radiation can alter the jet energies and multiplicities 
of the events, introducing distortions into the measured nf^°, 7?^“° and distributions. This 

effect is evaluated by performing pseudo-experiments using two dedicated signal samples generated with 
AcerMC [30] in combination with Pythia P2011C for hadronisation and parton showering. In these 
samples some Pythia P2011C parameters that control the showering are varied in ranges that are compat¬ 
ible with a study of additional jets in tt events [71], and half the difference of these two extremes is used 
as the systematic uncertainty. 


Underlying event and colour reconnection 

These systematic uncertainties are estimated using samples simulated with PowHEG-hvq and Pythia. The 
underlying-event uncertainty is obtained by comparing a sample with the Perugia 2012 tune (P2012) to a 
sample with the P2012 mpiHi tune [28]. The full difference in the fitted mass of the two models is taken 
as the systematic uncertainty for this source. Similarly, the colour reconnection systematic uncertainty is 
assigned as the difference in the fitted parameters of samples obtained with the P2012 and P2012 loCR 
tunes [28]. The same matrix-element-level PowHEG-hvq events generated with the CTIO PDFs are used 
for the three MC samples. The P2012 mpiHi tune is a variation of the P2012 tune with more semi- 
hard multiple parton interactions. The colour reconnection parameters were kept fixed to the P2012 tune 
values. Compared to the standard P2012 tune the P2012 loCR tune leads to significantly less activity in 
the transverse region with respect to the leading charged-particle as measured in Ref. [51]. In addition 
to assessing the effect of colour reconnection, this tune is therefore also used to estimate the systematic 
uncertainty associated with the particle spectra in the underlying event. 


Parton distribution functions 

The signal samples are generated using the CTIO PDFs. These PDFs, obtained from experimental data, 
have an uncertainty that is reflected in 26 pairs of possible PDF variations provided by the CTEQ group. 
To evaluate the impact of the PDF uncertainty on the tt signal templates, the events, from a sample gener¬ 
ated using MC@NFO with Herwig fragmentation, are re-weighted with the corresponding ratio of PDFs, 
and 26 pairs of signal templates are constructed, one pair per PDF uncertainty. For each pair, the average 
measured nitop is obtained from 500 pseudo-experiments each for the upward and downward variations 
of the PDF uncertainty. The corresponding uncertainty is obtained as half the difference of the two val¬ 
ues. From those the CTIO contribution is calculated as the sum in quadrature of the 26 uncertainties and 
amounts to 0.13 GeV and 0.10 GeV for the tt lepton-i-jets and tt dilepton analysis respectively. 

In addition, the signal tt samples are re-weighted to match the central PDFs for either the MSTW2008 [38] 
or the NNPDF23 [41] PDFs. The corresponding differences, taken as uncertainties, are 0.03 GeV and 
0.21 GeV for the tt lepton-i-jets analysis, and 0.01 GeV and 0.01 GeV for the tt dilepton ana¬ 
lysis. The final PDF systematic uncertainty is the sum in quadrature of the three contributions discussed 
above. 
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7.3. Modelling of non-tt processes 

The uncertainty in the modelling of non-tt processes is taken into account by varying the normalisation 
and the shape of the distributions of several contributions. 

The uncertainty on the Vk-tjets background determined from data [64] is dominated by the uncertainty on 
the heavy-flavour content of these events and amounts to ±30% of the overall normalisation. The same 
normalisation uncertainty is assigned to the Z±jets background normalisation. Uncertainties related to 
the VT±jets background shape are also considered. These stem from the variation of the heavy-flavour 
composition of the samples and from re-weightings of the distributions to match the predictions of Alp- 
gen. For the re-weighting, parameters are varied which affect the functional form of the factorisation 
and renormalisation scales, and the threshold for the matching scale used to connect the matrix-element 
calculation to the parton shower. 

The estimate of the background from NP/fake leptons determined from data is varied by ±50% to account 
for the uncertainty of this background source [65]. Uncertainties affecting the shape of this background 
are also included. For the NP/fake-electron background, the effects on the shape arising from the effi¬ 
ciency uncerfainfies for real and fake elecfrons are evaluafed and added in quadrafure. For fhe NP/fake- 
muon background, fwo differenl mafrix mefhods were used and averaged: fheir difference is faken as fhe 
sysfemafic uncerfainfy. 

In addition, the impact of changing the normalisation of the single top quark processes according to 
the uncertainty on the corresponding theoretical cross sections is considered. This yields a negligible 
systematic uncertainty in both the tt lepton±jets and tt —> dilepton analyses. 


7.4. Detector modelling 

Jet energy scale 

The JES is derived using information from test-beam data, LHC collision data, and simulation. The relat¬ 
ive JES uncertainty varies from about 1% to 3% depending on jet px and 77 as given in Ref. [58]. Since the 
estimation of the jet energy scale involves a number of steps, the JES uncertainty has various components 
originating from the calibration method, the calorimeter response, the detector simulation, and the specific 
choice of parameters in the physics model employed in the MC event generator. The total uncertainty is 
expressed in terms of 21 pj- and p-dependent components which are considered uncorrelated [58]. The 
uncertainties for the individual components and their sum are given in Table 4 in Appendix A. Despite the 
simultaneous fit of mtop, JSE and bJSE in the tt lepton±jets channel there is a non-negligible residual 
JES uncertainty. This is introduced by the variation of the jet energy scale corrections and their uncer¬ 
tainties with jet kinematics, which cannot be fully captured by global scale factors (JSE, bJSE). However 
the overall JES uncertainty is a factor of two smaller than in a one-dimensional analysis exploiting only 
templates of mj®™. In the tt dilepton channel, the contribution from the JES uncertainty constitutes the 
main component of systematic uncertainty on mtop. 

ft-jet energy scale 

This uncertainty is uncorrelated with the JES uncertainty and accounts for the remaining differences of 
f 7 -jefs and lighf-jefs affer fhe global JES was defermined. Eor fhis, an exfra uncertainly ranging from 
0.7% fo 1.8% and depending on jel px and p is assigned fo ft-jefs, due fo differences belween jefs con- 
laining fj-hadrons and fhe inclusive jel sample [58]. This additional sysfemafic uncerfainfy was obfained 
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from MC simulation and was verified using f?-tagged jets in data. The validation of the b-]&t energy 
scale uncertainty is based on the comparison of the jet transverse momentum as measured in the calor¬ 
imeter to the total transverse momentum of charged-particles associated with the jet. These transverse 
momenta are evaluated in the data and in MC simulated events for all jets and for fj-jets [58]. In addi¬ 
tion, a validation using tt lepton-f-jets events was performed. Effects stemming from f^-quark frag¬ 
mentation, hadronisation and underlying soft radiation were studied using different MC event generation 
models [58]. Thanks to the simultaneous fit to together with and the tt lepton-i-jets 
three-dimensional analysis method mitigates the impact of this uncertainty, and reduces it to 0.06 GeV, in¬ 
stead of 0.88 GeV in a two-dimensional analysis method (exploiting two-dimensional templates of 
and as in Ref. [8]), albeit at the cost of an additional statistical component of 0.67 GeV. In the 

tt —> dilepton channel, the contribution from the bJES uncertainty represents the second largest compon¬ 
ent of systematic uncertainty on mtop. 


Jet energy resolution 

To assess the impact of this uncertainty, before performing the event selection, the energy of each re¬ 
constructed jet in the simulation is smeared by a Gaussian function such that the width of the resulting 
Gaussian distribution corresponds to the one including the uncertainty on the jet energy resolution [72]. 
The fit is performed using smeared jets and the deviation from the central result is assigned as a systematic 
uncertainty. 


Jet reconstruction efficiency 

The jet reconstruction efficiency for data and the MC simulation is found to be in agreement with an 
accuracy of better than ±2% [73]. To account for the residual uncertainties, 2% of jets with < 30 GeV 
are randomly removed from MC simulated events. The event selection and the fit are repeated on the 
changed sample. The changes in the fitted parameters relative to the nominal MC sample are assigned as 
systematic uncertainty. 


Jet vertex fraction 

Residual differences between data and MC in the description of the fraction of the jet momentum asso¬ 
ciated with tracks from the primary vertex (used to suppress pile-up interactions) is corrected by apply¬ 
ing scale factors. These scale factors, varied according to their uncertainty, are applied to MC simula¬ 
tion events as a function of the jet pj. The resulting variation in the measured top quark mass in the 
tt —> lepton-i-jets analysis is 10 MeV, while it is negligible for the tt dilepton analysis. 

ft-tagging efficiency and mistag rate 

To account for potential mismodelling of the f?-tagging efficiency and the mistag rate, f?-tagging scale 
factors, together with their uncertainties, are derived per jet [61-63,74]. They are applied to the MC 
events and depend on the jet px and p and the underlying quark flavour. In this analysis these correction 
factors are obtained from dijet [62] and tt dilepton events. The same f7-tagging calibrations are applied 
to both the f^-i-jets and dilepton final states. The tt-based calibrations are obtained using the methodology 
described in Ref. [63], applied to the 7 TeV data. The statistical correlation stemming from the use 
of partially overlapping data sets for the tt dilepton mtop analysis and the fj-tagging calibration is 
estimated to be negligible. The correlation of those systematic uncertainties that are in common for the 
f?-tagging calibration and the present analyses is taken into account. Similarly to the JES uncertainty. 
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ff —> lepton+jets 

ff dilepton 

Combination 


KT [Gev] 

JSF 

bJSF 

[GeV] 

mfoT" [GeV] 

P 

Results 

172.33 

1.019 

1.003 

173.79 

172.99 


Statistics 

0.75 

0.003 

0.008 

0.54 

0.48 

0 

- Stat. comp, (mtop) 

0.23 

nja 

n/a 

0.54 



- Stat. comp. (JSF) 

0.25 

0.003 

n/a 

n/a 



- Stat. comp. (bJSF) 

0.67 

0.000 

0.008 

n/a 



Method 

0.11 +0.10 

0.001 

0.001 

0.09 + 0.07 

0.07 

0 

Signal MC 

0.22 + 0.21 

0.004 

0.002 

0.26 + 0.16 

0.24 

+ 1.00 

Hadronisation 

0.18 + 0.12 

0.007 

0.013 

0.53 + 0.09 

0.34 

+ 1.00 

ISR/FSR 

0.32 + 0.06 

0.017 

0.007 

0.47 + 0.05 

0.04 

-1.00 

Underlying event 

0.15 + 0.07 

0.001 

0.003 

0.05 + 0.05 

0.06 

-1.00 

Colour reconnection 

0.11 +0.07 

0.001 

0.002 

0.14 + 0.05 

0.01 

-1.00 

PDF 

0.25 + 0.00 

0.001 

0.002 

0.11+0.00 

0.17 

+0.57 

VT/Z+jets norm 

0.02 + 0.00 

0.000 

0.000 

0.01 + 0.00 

0.02 

+ 1.00 

W/Z+jets shape 

0.29 + 0.00 

0.000 

0.004 

0.00 + 0.00 

0.16 

0 

NP/fake-lepton norm. 

0.10 + 0.00 

0.000 

0.001 

0.04 + 0.00 

0.07 

+ 1.00 

NP/fake-lepton shape 

0.05 + 0.00 

0.000 

0.001 

0.01 + 0.00 

0.03 

+0.23 

Jet energy scale 

0.58 + 0.11 

0.018 

0.009 

0.75 + 0.08 

0.41 

-0.23 

/7-jet energy scale 

0.06 + 0.03 

0.000 

0.010 

0.68 + 0.02 

0.34 

+ 1.00 

Jet resolution 

0.22 + 0.11 

0.007 

0.001 

0.19 + 0.04 

0.03 

-1.00 

Jet efficiency 

0.12 + 0.00 

0.000 

0.002 

0.07 + 0.00 

0.10 

+ 1.00 

Jet vertex fraction 

0.01 + 0.00 

0.000 

0.000 

0.00 + 0.00 

0.00 

-1.00 

/7-tagging 

0.50 + 0.00 

0.001 

0.007 

0.07 + 0.00 

0.25 

-0.77 

^miss 

0.15 + 0.04 

0.000 

0.001 

0.04 + 0.03 

0.08 

-0.15 

Leptons 

0.04 + 0.00 

0.001 

0.001 

0.13 + 0.00 

0.05 

-0.34 

Pile-up 

0.02 + 0.01 

0.000 

0.000 

0.01 + 0.00 

0.01 

0 

Total 

1.27 + 0.33 

0.027 

0.024 

1.41 +0.24 

0.91 

-0.07 


Table 3: The measured values of mtop and the contributions of various sources to the uncertainty in the ti 
lepton+jets and the ff —> dilepton analyses. The corresponding uncertainties in the measured values of the JSF and 
bJSF are also shown for the ff —» lepton+jets analysis. The statistical uncertainties associated with these values 
are typically 0.001 or smaller. The result of the mtop combination is shown in the rightmost columns, together with 
the correlation (p) within each uncertainty group as described in Sect. 8. The symbol n/a stands for not applicable. 
Values quoted as 0.00 are smaller than 0.005. Finally, the last line refers to the sum in quadrature of the statistical 
and systematic uncertainty components. 


the uncertainty on the correction factors for the Zr-tagging efficiency is separated into ten uncorrelated 
components. The systematic uncertainty is assessed by changing the correction factor central values by 
±lcr for each component, and performing the fit. The final uncertainty due to the f?-tagging efficiency 
is calculated as the sum in quadrature of all contributions. A similar procedure is applied for the mistag 
rates for c-jets, albeit using four separate components. In addition, the correction factors and mistag rates 
for light-jets are varied within their uncertainty, and the corresponding shifts in the measured quantities 
are summed in quadrature. The size of the fi-tagging systematic uncertainty of 0.50 GeV observed in the 
tt —> lepton+jets analysis is mostly driven by the induced change in shape of the distribution. 

Lepton momentum and missing transverse momentum 

The lepton momentum and the are used in the event selection and reconstruction. For the leptons, the 
momentum scale, resolution and identification efficiency are measured using high-purity Z ^ ii data [54, 
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55]. The uncertainty due to any possible miscalibration is propagated to the analyses by changing the 
measured reconstruction efficiency, lepton pj, and the corresponding resolution, within uncertainties. 

The uncertainties from the energy scale and resolution corrections for leptons and jets are propagated 
to the S™®®. The systematic uncertainty related to the accounts for uncertainties in the energies 
of calorimeter cells not associated with the reconstructed objects, and from cells associated with low-pT 
jets (7 GeV< pj < 20 GeV), as well as for the dependence of their energy on the number of pile-up 
interactions [60]. 


Pile-up 

The residual systematic uncertainty due to pile-up was assessed by determining the dependence of the 
fitted top quark mass on the amount of pile-up activity, combined with uncertainties in modelling the 
amount of pile-up in the sample. 


7.5. Summary 

The resulting sizes of all uncertainties and their sum in quadrature are given in Table 3. The total un¬ 
certainties on JSF, bJSF and amount to 1.27 GeV, 0.027, 0.024 and 1.41 GeV, respectively. 

Within uncertainties, the fitted values of JSF and bJSF are consistent with unity. 


8. Combination of the mtop results 


The results of the tt —> lepton-i-jets and tl —> dilepton analyses listed in Table 3 are combined using 
the Best Linear Unbiased Estimate (BLUE) method [75,76], implemented as described in Refs. [77, 
78]. The BLUE method determines the coefficients (weights) to be used in a linear combination of the 
input measurements by minimising the total uncertainty of the combined result. In the algorithm, both 
the statistical and systematic uncertainties, and the correlations (p) of the measurements, are taken into 
account, while assuming that all uncertainties are distributed according to Gaussian probability density 
functions. 


8.1. Correlation of the tt lepton-i-jets and tt dilepton measurements 

To perform the combination, for each source of systematic uncertainty, the uncertainties as well as the 
correlation of the measurements of mtop were evaluated. 

The measurements are taken as uncorrelated for the statistical, the method calibration and the pile-up 
uncertainties. Eor the remaining uncertainty components there are two possible situations. Either the 
measurements are fully correlated, p = +1, i.e. a simultaneous upward variation of the systematic un¬ 
certainty results in a positive (or negative) shift of mtop for both measurements, or fully anti-correlated, 
p - - \.\n the latter case one measurement exhibits a positive shift and the other a negative one. 

Eigure 6(a) shows the two dimensional distribution of the systematic uncertainties, denoted by 

and obtained in the ^-i-jets and dilepton analyses for all components of the sources of systematic 
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Figure 6; The systematic uncertainties of mtop in the f+jets analysis versus those of the dilepton analysis. Fig¬ 
ures (a-c) refer to the results evaluated for the three-dimensional analysis (3d), two-dimensional analysis (2d) and 
one-dimensional analysis (Id). The points show the estimated systematic uncertainties on mtop for the two analyses, 
and the uncertainty bars reflect the corresponding statistical uncertainties. The different colours reflect the different 
correlations described in Sect. 8.1. 


uncertainty for which the measurements are correlated. The points show the estimated size of the un¬ 
certainties, and the error bars represent the statistical uncertainties on the estimates. Some uncertainty 
sources in Table 3, such as the uncertainty related to the choice of MC generator for signal events, con¬ 
tain only a single component. For these type of sources, the correlation is either p - +\ (red points) or 
p - (blue points). The size of the uncertainty bars in Fig. 6(a) indicates that the distinction between 
p - +\ and p = -1 can be unambiguously made for all components that significantly contribute to the 
systematic uncertainty on mtop. For uncertainty sources that contain multiple components such as the JES 
uncertainty described in Appendix A, the correlations given in Table 3 differ from p = +1. For fhese cases 
fhe correlation is obfained by adding fhe corresponding covariance terms of fhe componenfs and dividing 
by fhe respective fofal uncerfainfies of fhe source. 

For each sysfemafic uncerfainfy, fhe size of and Am^^p, and fhe correlafion of fhe measuremenfs 

depend on fhe defails of fhe analyses. This can be seen from Figs. 6(b, c) where fhe same informafion as 
in Fig. 6(a) is shown, buf for differenl implemenfafions of fhe /’-i-jels analysis, while leaving fhe dilepfon 
analysis unchanged. Figure 6(b) corresponds fo a fwo-dimensional analysis, similar fo Ref. [8], which is 
realised by fixing fhe bJSF fo unify. Finally, Fig. 6(c) shows fhe resulf of a one-dimensional analysis, in 
which fhe values of fhe JSF and bJSF are fixed fo unify. For fhis implemenfafion, as for fhe dilepfon 
analysis, only mtop is obfained from fhe fif fo dafa. Compared fo fhe fwo-dimensional analysis, fhe 
fhree-dimensional analysis reduces some sources of uncertainly on m^p. As an example, fhe righfmosf 
red poinf in Fig. 6(b), which corresponds fo fhe bJES uncerfainfy, lies close fo fhe vertical line in Eig. 6(a), 
i.e. for fhe ^-i-jels analysis fhe impacf of fhis source was considerably reduced by fhe bJSE determination 
from dafa. The change in fhe correlalions of fhe measuremenfs for specific sources of uncerfainfy, caused 
by a variation of fhe analysis sfrafegy, is apparenf from Eig. 6(c), where for bofh analyses only m^p is 
obfained from fhe dafa. In fhis case fhe exploiled observables are much more similar and consequenlly, 
fhe measuremenfs of mtop are fully correlated for all sources of uncerfainfy lhal significanlly conlrib- 
ule fo fhe fofal uncertainly. This demonslrafes lhal fhe fhree-dimensional analysis nol only reduces fhe 
impacf of some sources of uncertainly, mainly fhe JES and bJES uncerfainfies, buf also makes fhe Iwo 
measuremenfs less correlaled, Ihus increasing fhe gain in fhe combination of fhe fwo esfimales of mtop. 

To besl profil from fhe combinalion of fhe fwo measuremenfs, Iheir correlafion should be as small as 
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possible, see Ref. [78]. Consequently, the jet energy scale factors measured in the f’+jets analysis have 
not been propagated to the dilepton analysis, as was first done in Ref. [79]. Transferring the scales would 
require adding an additional systematic uncertainty to the dilepton analysis to account for the different 
jet energy scale factors caused by different kinematical selections and jet topologies of the two analyses. 
The two final sfafes confain eifher fwo or four jefs fhaf have differenf disfribufions in jef pT> and differ- 
enf amounfs of final sfafe QCD radiation. Mosf nofably, fhis would also resulf in a large correlation of 
fhe measuremenfs, similar fo fhaf observed for fhe one-dimensional analyses shown in Fig. 6(c). Con- 
sequenfly, fhe knowledge of nitop from fhe T’-i-jefs analysis would nol significanfly improve when including 
a dilepfon measuremenf obfained wifh Iransferred jef energy scales. For an example of such a sifuafion 
see Table VI of Ref. [79]. 

Using fhe correlafions defermined above, fhe combinafion of fhe mtop resulfs of fhe tt —> lepfon-i-jefs and 
tt —> dilepfon analyses yields: 


= 172.99 ± 0.48 (sfaf) ± 0.78 (sysf) GeV = 172.99 ± 0.91 GeV. 

This value corresponds fo a 28% gain in precision wifh respecf fo fhe more precise T’-f-jefs measuremenf. 
The compafibilify of fhe inpuf measuremenfs is very good, and corresponds fo 0.750" = 

-1.47 ± 1.96 GeV). The BLUE weighfs of fhe resulfs of fhe tt lepfon-i-jefs and tt —> dilepfon analyses 
are 54.8% and 45.2%, respectively. The fofal correlation of fhe inpuf measuremenfs is -7% and fhe;y^ 
probabilify of fhe combinafion is 45.5%. The lisf of all uncerfainfies of fhe combined resulf, fogefher wifh 
fhe correlation of fhe measuremenfs for each group of uncerfainfies, is provided in Table 3. The currenf 
precision is mosfly limited by sysfemafic uncerfainfies relafed fo fhe MC modelling of tt evenfs, and fo 
fhe calibration of fhe jef energy scales. 


8.2. Stability of the results 

The dependence of fhe combined resulf on fhe sfafisfical uncerfainfies of fhe evaluated sysfemafic un¬ 
cerfainfies is invesfigafed by performing one fhousand BLUE combinafions in which all inpuf uncer¬ 
fainfies are independenfly smeared using Gaussian functions cenfred af fhe expecfed values, and wifh a 
widfh corresponding fo fheir sfafisfical uncerfainfies. Using fhe smeared uncerfainfies, fhe correlafions 
are re-evaluafed for each pseudo-experimenf. The combined nitop and ifs fofal uncerfainfy are disfrib- 
ufed according fo Gaussian funclions of widfh 37 MeV and 43 MeV, respectively. Similarly, fhe BLUE 
combinafion weighfs and fhe fofal correlation are Gaussian disfribufed, wifh widfhs of 2.5% and 6.1%, 
respecfively. These effecls are found fo be negligible compared fo fhe fofal uncertainly of fhe combined 
resulf. Consequenlly, no addifional sysfemafic uncertainly is assigned. 


9. Conclusion 

The lop quark mass was measured via a fhree-dimensional femplale melhod in fhe tt leplon-i-jels final 
sfafe, and using a one-dimensional femplale melhod in fhe ti dilepfon channel. Bolh analyses are 
based on yfs = 7 TeV prolon-prolon collision ATLAS dafa from fhe 2011 LHC run corresponding fo 
an integrated luminosity of 4.6 fb“^ In fhe f’-f-jels analysis, mtop is determined fogefher wifh a global jef 
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energy scale factor (JSF) and a residual Zr-to-light-jet energy scale factor (bJSF). The measured values 
are: 


= 172.33 ± 0.75 (stat + JSF + bJSF) ± 1.02 (syst) GeV, 
JSF = 1.019 ± 0.003 (stat) ± 0.027 (syst), 
bJSF = 1.003 + 0.008 (stat) ± 0.023 (syst), 

= 173.79 ± 0.54 (stat) ± 1.30 (syst) GeV. 


These measurements are consistent with the ATLAS measurement in the fully hadronic decay chan¬ 
nel [13], and supersede the previous result described in Ref. [8]. 

A combination of the tt lepton-i-jets and tt —> dilepton results is performed using the BLUE technique, 
exploiting the full uncertainty breakdown, and taking into account the correlation of the measurements 
for all sources of the systematic uncertainty. The result is: 


= 172.99 ± 0.48 (stat) ± 0.78 (syst) GeV = 172.99 ± 0.91 GeV. 

This corresponds to a gain in precision with respect to the more precise ^-f-jets measurement of 28%. The 
total uncertainty of the combination corresponds to 0.91 GeV and is currently dominated by systematic 
uncertainties due to jet calibration and modelling of the tt events. 
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Appendix 

A. Jet energy scale uncertainty: detailed components 

The relative JES uncertainty varies from aboutabout 1% to 3% depending on jet properties as given in 
Section 13 of Ref. [58]. These components correspond to the eigenvectors of the reduced covariance 
matrix for the JES uncertainties, as described in Section 13.3 of Ref. [58]. The initial sources of nuisance 
parameters (NP) originating from the in-situ determination of the JES are listed in Table 10 of Ref. [58]. 
According to their nature, they are categorised into the classes: detector description, physics modelling, 
statistics and method, mixed detector and modelling. Einally, following Section 13.6 of Ref. [58], a re¬ 
duction of the number of nuisance parameters is performed for each category giving various components. 
Their px dependences are given in Eig. 42 of Ref. [58]. The total JES uncertainty is provided together 
with its 21 sub-components in Table 4. Their separate effects on the fitted top quark mass are summed in 
quadrature to determine the total jet energy scale uncertainty given in Table 3. Eor further details about 
each component, see Ref. [58]. 
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ti —> lepton+jets 

tt —> dilepton 

Combination 


[GeV] 

AJSF 

AbJSF 

Amf“ [GeV] 

Au^^^iGeV] 

P 

Statistical (total) 

0.18 ±0.04 

0.003 

0.001 

0.16 ±0.03 

0.11 

-0.25 

- Statistical NPl 

-0.17 ±0.02 

±0.002 

±0.001 

±0.01 ± 0.02 

0.09 

-1.00 

- Statistical NP2 

±0.02 ± 0.00 

±0.001 

-0.000 

±0.05 ± 0.00 

0.03 

±1.00 

- Statistical NP3 

-0.01 ± 0.02 

±0.001 

±0.001 

±0.12 ±0.02 

0.05 

-1.00 

- rj inter-calibration (stat.) 

-0.07 ± 0.02 

±0.001 

±0.001 

±0.10 ± 0.02 

0.01 

-1.00 

Modelling (total) 

0.31 ±0.06 

0.009 

0.002 

0.52 ± 0.04 

0.26 

-0.18 

-Modelling NPl 

-0.30 ± 0.03 

±0.006 

±0.001 

±0.22 ± 0.02 

0.07 

-1.00 

- Modelling NP2 

±0.03 ± 0.02 

±0.002 

-0.000 

±0.14 ±0.02 

0.08 

±1.00 

- Modelling NP3 

-0.01 ± 0.02 

-0.002 

-0.000 

-0.15 ±0.02 

0.07 

±1.00 

- Modelling NP4 

-0.01 ± 0.00 

±0.000 

±0.000 

±0.02 ± 0.00 

0.00 

-1.00 

- T) inter-calibration (model) 

±0.07 ± 0.04 

±0.007 

-0.001 

±0.43 ± 0.03 

0.23 

±1.00 

Detector (total) 

0.05 ± 0.03 

0.007 

0.001 

0.45 ± 0.04 

0.20 

-0.19 

- Detector NPl 

-0.01 ±0.03 

±0.007 

±0.001 

±0.45 ± 0.02 

0.20 

-1.00 

- Detector NP2 

-0.05 ± 0.00 

±0.000 

±0.001 

±0.03 ± 0.00 

0.02 

-1.00 

Mixed (total) 

0.02 ± 0.02 

0.001 

0.001 

0.03 ± 0.02 

0.01 

-0.80 

-Mixed NPl 

-0.02 ± 0.00 

±0.000 

±0.001 

±0.02 ± 0.00 

0.00 

-1.00 

- Mixed NP2 

±0.00 ± 0.02 

±0.001 

-0.000 

±0.02 ± 0.02 

0.01 

±1.00 

Single particle high-pT 

±0.00 ± 0.00 

±0.000 

-0.000 

±0.00 ± 0.00 

0.00 

±1.00 

Relative non-closure MC 

±0.00 ± 0.02 

±0.001 

-0.000 

±0.03 ± 0.02 

0.02 

±1.00 

Pile-up (total) 

0.15 ±0.04 

0.001 

0.002 

0.04 ± 0.03 

0.09 

±0.03 

-Pile-up: Offset(p) 

-0.11 ±0.02 

-0.001 

±0.001 

-0.02 ± 0.02 

0.07 

±1.00 

-Pile-up: Offset(nvtx) 

-0.10 ±0.04 

-0.000 

±0.001 

±0.03 ± 0.03 

0.04 

-1.00 

Flavour (total) 

0.36 ± 0.04 

0.012 

0.008 

0.03 ± 0.03 

0.20 

-0.17 

- Flavour Composition 

-0.24 ± 0.02 

±0.006 

-0.002 

-0.02 ± 0.02 

0.14 

±1.00 

- Flavour Response 

-0.28 ± 0.03 

±0.011 

-0.008 

±0.03 ± 0.02 

0.14 

-1.00 

Close-by jets 

-0.22 ± 0.04 

±0.005 

±0.002 

±0.25 ± 0.03 

0.01 

-1.00 

b-ieX energy scale 

±0.06 ± 0.03 

±0.000 

±0.010 

±0.68 ± 0.02 

0.34 

±1.00 

Total (without bJES) 

0.58 ±0.11 

0.018 

0.009 

0.75 ± 0.08 

0.41 

-0.23 


Table 4: The individual components of the JES uncertainty according to Ref. [58], together with the corresponding 
uncertainties on JSF, bJSF, and Some components listed are calculated as the sum in quadrature 

of several sub-components. The corresponding measurement correlations per group described in Sect. 8 are also 
reported. 
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